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Abstract 



We outline a method based on successive canonical transformations which yields 
a product expansion for the evolution operator of a general (possibly non-Hermitian) 
Hamiltonian. For a class of such Hamiltonians this expansion involves a finite num- 
ber of terms, and our method gives the exact solution of the corresponding time- 
dependent Schrodinger equation. We apply this method to study the dynamics 
of a general nondegenerate two-level quantum system, a time-dependent classical 
harmonic oscillator, and a degenerate system consisting of a spin 1 particle in- 
teracting with a time-dependent electric field £ (t) through the Stark Hamiltonian 



Recently, a method based on successive canonical transformations has been used to obtain 
exact solution of the Schrodinger equation 



for a class of dipole Hamiltonians [T], 0, [| and time-dependent harmonic oscillators 0. 
For these systems the Hamiltonian is a nondegenrate Hermitian operator. The purpose 
of the present article is to extend the application of this method to the cases where the 
Hamiltonian is non-Hermitian and involves degenerate eigenvalues. 

Non-Hermitian Hamiltonians have been used to model a variety of physical systems in- 
volving decaying states, ||. The solution of the Schrodinger equation for a time-dependent 
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H = X(J-£) 2 . 



I. Introduction 



i-m)) = mm*)) 




two- level non-Hermitian Hamiltonian has been considered in Refs. |6], [/J. Another mo- 
tivations for the study of the Schrodinger equation for a time-dependent non-Hermitian 
Hamiltonian is the fact that the solution of every linear ordinary differential equation 
(ODE) may be reduced to the solution of a system of first order linear ODEs which can 
be written in the form of the time-dependent Schrodinger equation (|l|) or alternatively 



where U(t) is the evolution operator. For a general linear ODE the corresponding Hamil- 
tonian H(t) may be a non-Hermitian matrix with degenerate eigenvalues. 

The method of adiabatic product expansion developed in Refs. H does not di- 
rectly apply to quantum systems with non-Hermitian Hamiltonians. In this article we 
shall present a generalization of this method which applies to arbitrary (possibly) non- 
Hermitian Hamiltonians with degenerate as well as nondegenerate eigenvalues. 

The organization of the article is as follows. In section II we review the basic results 
concerning the adiabatic approximation for degenerate and non-Hermitian Hamiltonians. 
In section III we discuss the generalization of the method of adiabatic product expansion 
to these Hamiltonians. In section IV we use the results of section III to study the so- 
lution of the Schrodinger equation for a general nondegenerate non-Hermitian two-level 
Hamiltonian. In section V, we apply the general results obtained in section IV to treat 
the classical equation of motion for a harmonic oscillator with a time-dependent fre- 
quency. In section VI, we discuss the application of the adiabatic product expansion to 
study the quadrupole interaction of a spin 1 particle with a time- dependent electric field 
£ = (£i(t), £2(t), 0). We show that the corresponding Hamiltonian which has a degenerate 
and a nondegenerate eigenvalue is canonically equivalent to a Hamiltonian which has only 
nondegenerate eigenvalues. Furthermore, we show that if the direction of the electric field 
depends in a particular way on its magnitude, then our method yields the exact solution 
of the Schrodinger equation. Finally we present our conclusions in section VII. 

II. Adiabatic Approximation for Non-Hermitian Hamil- 
tonians 

Let H = H[R] be a parametric Hamiltonian which depends on a set of real parameters 
R = (R 1 , R 2 , ■ • • , R d ) labelling the points of a smooth manifold M. Let E n [R] denote the 
eigenvalues of H[R] and TC n [R] be the degeneracy subspace associated with E n [R]. Let A/" 
denote the degree of degeneracy of E n [R], i.e., the complex dimension of H n [R). We shall 
assume that the spectrum of H [R] is discrete and A/" does not depend on R. 
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Now let \ip n ,a;R) and \<f) n ,a;R) form a complete biorthonormal basis of the Hilbert 
space H This means that |-0 n , a; i?) with a G {1, 2, • ■ ■ ,7V} form a basis of H n [R], in 
particular 

fl-[i2]|V„,o;i2) = S n [i2]|^,o;i2) , (5) 

and \<p n , a;R) satisfy 

H[R]*\<j> n ,a-,R) = E*[R]\<j> n ,a-,R) , (6) 

(4> m , b; R\ip n , a; R) = 8 mn 5 ab , (7) 
W 

EE> n ,a;i2)<0 n ,a;i2| = l. (8) 

n a=l 

Next suppose that the parameters R % depend on time t, then R(t) defines a curve 
C in the parameter space M, and the Hamiltonian, its eigenvalues and eigenvectors 
become time-dependent. In this case we use the notation H(t) := H[R(t)}, E„(t) : = 
E n [R(t)}, \ip n ,a;t) := \ip n , a; R(t)) , and \<f) n ,a;t) := \<f) n ,a;R(t)). We shall assume that 
E n (t), \i/) n ,a;t) and \<p n , a;t) are smooth functions of t and that during the evolution of 
the system the eigenvalues of the Hamiltonian do not cross, i.e., if E m (0) < E n (0), then 
for all t G [0,t], E m {t) < E n (t), where r denotes the duration of the evolution of the 
system. 

Differentiating both sides of Eq. ([5]) with respect to t, taking the inner product of both 
sides of the resulting equation with \<f> m , b; t), for arbitrary m and b, and using Eqs. - 
(0), we have 

\E m (t) - E n (t)](<f> m , b; A-\^ n: a; t) + (<j> m , b; t\H(t)\il> n , a; t) - 5 mn 5 ab E(t) = . (9) 
Here a dot denotes differentiation with respect to t. For m ^ n, Eq. @ reads 

Now let us express the solution of the Schrodinger equation (P in the basis {\il> n , a; t)}. 
Then 

TV 

l#)} = EECWIi,M), (ii) 

n a=l 

where C™(t) are complex coefficients. Substituting Eq. ( |TTD in the Schrodinger equa- 
tion (HD, taking the inner product of both sides of the resulting equation with \<j) m ,b;t), 
and making use of Eqs. (||), @, (|7|), and fllOD , we find 

iCr - E m C™ + E UK, b; t\^ m , a; t)C? = -i £ £ . (12) 

o=l dt n^ma=l E n (t) - E m (t) 

The special case of this equation with Af = 1, i.e., the nondegenerate case, has been origi- 
nally derived by Garrison and Wright |J in their investigation of the adiabatic geometric 
phase ]TD for non-Hermitian Hamiltonians H |TT, I2[ 13 ]. 



If the right-hand side of Eq. fll2"|) is negligible, then one says that the system undergoes 
an adiabatic evolution, (15], [16], Q, [L7|] . In this case, the equations for C™ decouple and 
their solution is given by 

TV 

QW = EW"(o)» (13) 

6=1 

where K™ b (t) are entries of the invertible matrix 

K n {t) :=e-^o En{s)ds Vexp t f m A n [R] , (14) 
V denotes the path-ordering operator, A n is the matrix of one-forms with entries 

A n ab [R] :=i(<p n ,a;R\d\i(> n ,b;R) , (15) 

d stands for the exterior derivative with respect to R\ and the line intergral in Eq. ( |T4"P is 
evaluated along the curve C defined by R(t). If C is a closed curve in M, the Hamiltonian 
has a periodic time- dependence and the path-ordered exponential in Eq. (|1"4"D, which takes 
the form 



Pexp 



(16) 



is the non-Hermitian analogue of the non-Abelian adiabatic geometric phase (18 



Note that if the initial vector \ip(0)) is an eigenvector of the initial Hamiltonian H(0), 
then the adiabaticity of the evolution implies that \ip(t)) is an eigenvector of H(t) for all 
t G [0, t}. In terms of the time-evolution operator U(t) of Eq. (|3]) this is expressed by 



U(t) 
U^(t) 



U (0) (t) , where 
N 

E E K b {t)\^a-t){^b-M 



(17) 
(18) 



n a, 6=1 



One can easily show that U^(t) is invertible, and its inverse is given by 



A'' 



^ (0) "W =E E K^\t)\^a^){(j> n ,b-t\ 



(19) 



n a, 6=1 

where K n 1 (t) is the inverse of K n (t). 



III. Adiabatic Canonical Transformations and the Gen- 
eralized Adiabatic Product Expansion 

Let g(t) be an invertible linear operator acting on the Hilbert space. Then the transfor- 
mations: 



\m) - Wit)):=g{t)m)) , 
H(t) - H\t) :=g(t)H(t)g(t)- 1 - 
U(t) -> U'(t):=g(t)U(t)g(0)- 1 , 



(20) 
(21) 
(22) 



leave the form of the Schrodinger equation invariant. We shall call such a transformation 
a canonical transformation. 

Now let us investigate the consequences of the canonical transformation defined by 
g(t) = £A°)(t) _1 . We shall call this transformation the adiabatic canonical transformation. 
Denoting the transformed Hamiltonian H' by H W ? we have 

TV M 

H {1) (t) = ]T Y,T, H ^ nm ( t Mn,a;0)(<f )m ,b ] 0\ , where (23) 

n,m^n 0=1 6=1 

#£ )nm (*) := -K c {ty l A™{t)K2{t) , and ^(f) := ifa, c; t\^ m , d; t) .(24) 

Because g(0) = C/^(0) _1 = 1, the transformed evolution operator is given by 

U'(t) = U i0) (t)- 1 U(t) . (25) 

Clearly if the adiabatic approximation is valid, H^'(t) ~ and U'(t) ~ 1. 

Let us suppose that H^ l > (t) has a discrete spectrum and denote by E$ (t) and Mi the 
eigenvalues of H^(t) and their degree of degeneracy. Furthermore, let {\ip!$, ai',t), \(j>£\ ai',t}} 
be a biorthonormal eigenbasis of the Hilbert space, i.e., 

H^(t)\^, ai ;t) = E^[R]\^, ai ;t) , 

EE ltfUi;i)(ft«i;t| = i. 

ni ai=l 

Then H^(t) shares the properties of the original Hamiltonian H(t), and we can repeat 
the above analysis using H^(t) in place of Hit). In this way the adiabatic approximation 
yields the approximate evolution operator 

Mi 

U {l \t) = Y. E K^i^mSl^^nlM-A (26) 

ni Oi,6i=l 

for H®(t), where K^(t) are the entries of the matrix K^ ni obtained by replacing 
E n (t), \ij n ,a;t) and \<f> n ,a;t) in Eqs. (H) and (0) by 1^, oi; t), and |0g>, au*), 

respectively. 

Next we perform the adiabatic canonical transformation defined by g(t) = tA 1 )^) -1 . 
This leads to a transformed Hamiltonian H ^ (£) which is related to fA 1 -* (t) according to 
Eqs. © and ([§) with if", |V>n,a;t) and \<j> n ,b',t) replaced by K^" 1 , \ip^,a x \t) and 
\4>^, a>i',t). The transformed evolution operator is given by 

u^ity^ity'uit) . 



Repeating this procedure we obtain, after N successive adiabatic canonical transfor- 
mations, a transformed Hamiltonian (£) and a transformed evolution operator which 
is given by 

^-1)^-1^-2)^-1 . _ _ uWty-iuQ . 

Here U^(t), with £ 6 {1, 2, • • • , N — 1}, denotes the approximate evolution operator 
obtained by performing adiabatic approximation on the Hamiltonian H^\t). 

If for some N the adiabatic approximation yields the exact solution of the Schrodinger 
equation for the Hamiltonian (£), then by construction H^ N+1 ^ (t) = and f/( Ar+1 ) (t) = 
1. In this case, the original evolution operator is given by 

U(t) = (fill® (*)••• U {N) (t) . (27) 

If the adiabatic approximation fails for all H^ N \t), then there are two possibilities: 

i) one obtains an infinite product expansion for the evolution operator 

oo 

U(t) = l[U (£) (t) :=U (0) (t)U (1 \t)---U (£) (t)-- - . (28) 

In this case, one may view Eq. (p7|) as a generalization of the adiabatic approxima- 
tion. 

ii) one obtains H®(t) = H^\t) for some i and j with i ^ j. In this case a direct 
application of the method of adiabatic product expansion does not produce a so- 
lution. However, as we shall see in the following section, sometimes it is possible 
to modify this method by combining the adiabatic canonical transformation with 
other canonical transformations, so that one obtains a finite or an infinite product 
expansion with distinct terms. 



IV. Application to Two-Level Hamiltonians 

Two-level nondegenerate Hamiltonians provide the simplest nontrivial quantum systems. 
This has been one of the main reasons for the study of these Hamiltonians since the 
early days of quantum mechanics. In this section we shall consider the most general 
nondegenerate two-level Hamiltonian which may or may not be Hermitian. 

In an arbitrary basis of the Hilbert space (C 2 ), the Hamiltonian is given by a two- by- 
two complex matrix H. One can perform a quantum canonical transformation (|2l|) defined 
by g{t) = exp{i /o[tr H(s)]ds/2} to map the Hamiltonian H to a traceless Hamiltonian 
of the form 

, a b , 

H:= , (29) 



where tr H denotes the trace of H, and a = a(t), b = b(t), c = c(t) are complex- valued 
smooth functions of t. 

We can easily solve the eigenvalue problem for the Hamiltonian (p9j). The eigenvalues 
are given by 



E x (t) := -E(t), E 2 (t):=E(t), where E := Va 2 + be. 



(30) 



We shall demand that during the time interval [0, r] of interest E 7^ 0, so that the 
eigenvalues are nondegenerate. In particular, no level crossings occur. Then a possible 
choice for a biorthonormal eigenbasis is 



-b 

a + E 



a + E 
c 



m ' ' N*[ a * + E* ) ' ' N* [ b* 

where R = (a, 6, c) and N := 2E(a + £)■ 

Next we compute £/ {0) and Using Eqs. (|1SD , (HD, (ID, and (||), we find 

= e(t)^i;O><^;O| + C(t)|V2;O>(0i;O|, 



(31) 
(32) 



(33) 
(34) 



where 



K\t) 
K\t) 

T](t) 

m 

a{t) 



K 2 u (t)=exp 



2 f E(s)ds 
Jo 

rl2 



'ir](t) /•*(*) 

2 iii(o) 

'— i?7(t) /*(*) 

2 7j?(o) 



(2E)- 1 {da + dE + 



cdb , 



(2E)- i (da + rfE + 



bde 
a + E 



Him 



mm 



ie 



-2ia(t) 



2 

^g2io(t) 



1 + 



a(t) 

W) 



d 
It 



c(t) 



1 + 



a(t) 
W) 



d_ 
~dt 



a(t) + E(t) 
b(t) 



a(t) + E(t) 



r)(t) i f R (V cdb — bde 



2 AJr( ) E(E + a) 
The transformed Hamiltonian has the following matrix expression 



(35) 
(36) 

(37) 
(38) 



HM(t) 



c (i)(t) - a «(t) 



(39) 



where 



a«(t) 



&o£(*) + CpC(t) 
2£ n 



(40) 



6(1)(i) ^ _ M«y|±|m , (41) 

m m _ -(a + ^o) 2 e(t) + c^(t) , . 

lJ 2£ (£ + a ) ' ^ 

a := a(O), 6 := 6(0) , c := c(0) , £ := £(0) , 

and we have used Eqs. fl3l"D , (]32f) , and 

Note that the transformed Hamiltonian if^ 1 -* (i) is traceless, and one can obtain H^{t) 



by substituting for a, for 6, for c, and -E^ 1 -* := y (aW) 2 + feWc^ 1 ) for in 
Eqs. (pip, and (^) - ([38]). Clearly this can be repeated indefinitely, and one can compute 
H(Q for arbitrary I. 

The adiabatic approximation corresponds to the cases where the matrix elements of 
can be neglected. As seen from Eqs. (|39p - (|42|) this happens whenever both £ 
and ( are negligible. One can also check that if only one of these quantities is negligible, 
then H^(t) is equal to the other times a constant matrix. This means that HW(t) 
has essentially stationary eigenvectors and the adiabatic approximation would yield the 
solution of the Schrodinger equation for H^ 1 ' . In fact, it is not difficult to check that for 
the cases that either £ or £ is negligible, H^ 2 \t) m 0. In particular, setting £ = or £ = 
implies (t) = and the evolution operator is given by 

U(t) = UW(t)UW(t) . (43) 

Therefore, the conditions £ = and £ = each define a class of exactly solvable two-level 
systems. In view of Eqs. (|36| ) and (|37D , these are 

Class 1: The two level systems for which = \x = constant, or alternatively 
c = fi(fib + y/4a 2 + n' 2 b 2 ) /2; 

Class 2: The two level systems for which — v — constant, or alternatively 
c = b/u 2 — a 2 /b. 

In general £ and £ do not vanish and the adiabatic product expansion does not termi- 
nate. There is also a special class of two-level systems for which the product expansion 
has a periodic structure in the sense of case (ii) of the preceding section. This is 

Class 3: The two level systems for which a = 0. 

Setting a = in Eqs. @, and (0) and defining /(£):= iy/c{t)/b{t), we have 

fof(t)e- iv ® H ^ _ /(*)e*rt*> 



where fo := /(O). Substituting these equations in Eq. (j3~9| ), we obtain 

= *■>(«) f ;° S " (t ' ) , where = ^ . (44) 

y f smr](t) - cos r)(t) J 4/W 

This Hamiltonian has two interesting properties. 
1. If 60 = Co, then fo = i and 

= £ (1) (t) [sin77(t)a 2 + cos77(t)<7 3 ] = ^(ty^/Vae"^ 1 / 2 , (45) 
where cij are Pauli matrices, and we have used the identity 

3 

e-^aje^ = co8{2(p)a j + sin(2y?) ^ijk^k , for i ^ k . (46) 

k=l 

In Eq. (Pp, </? is an arbitrary complex variable and is the totally anti-symmetric 
Levi Civita symbol with £123 = 1. For the time periods during which \Jb(t)c(t) 
is real, 77 (t) is real, and the Hamiltonian (fl"5| ) is anti-Hermitian. In particular, 
its eigenvectors are orthogonal. Up to a factor of i this Hamiltonian describes 
the interaction of a spin 1/2 magnetic dipole with a changing magnetic field. This 
system has a SU(2) dynamical group |L9], [j], H- For the time periods during which 



b(t)c(t) is imaginary, 77 (£) is imaginary, and up to a factor of i the Hamiltonian 
3D describes a quantum system with a ^[/(l, 1) dynamical group. A Hermitian 
analogue of such a system is the time-dependent generalized harmonic oscillator 



2. Performing the adiabatic canonical transformation on ([45"D , we arrive at the unex- 
pected result 

H^\t) = H{t). (47) 

Therefore, direct application of the method of adiabatic product expansion does not 
lead to a solution. 

Next we shall describe a modification of the method of adiabatic product expansion 
which yields an infinite product expansion for the evolution operator of the Class 3 systems 
which involve distinct terms. 

Consider the transformed Hamiltonian We can express this Hamiltonian using 

Eq. (|3UD with 

a (1) (t) = (t) cos V (t) , (i) = /"^W (t) sin 77(f) , c (1) (t) = / £ (1) (*) sin V (t) . 

(48) 

Although this Hamiltonian does not belong to Class 3, it can be canonically transformed 
to a Hamiltonian which belongs to Class 3, i.e., its diagonal matrix elements vanish. This 



transformation is defined by g(t) = exp{i J a^(s)dsas}. The corresponding transformed 
Hamiltonian is given by 

&i(t) 



Hi(t) 



ci(t) 



(49) 



where 



bJt) := b (1 \t)e^ {t) , ci(t) := c«(f)e-^ , and 7l (f) := 2 / a«(s)ds . (50) 

JO 

The evolution operator ?7i of f/i is related to the evolution operator of the original Hamil- 
tonian H according to 



Ui(t) = e^o^^ 3 ^ U (0 \ty U(t) 



(51) 



where we have used Eq. ([2"2]). 

Now since Hi has the same form as H, we can repeat the above analysis using H\ 
in place of H. Performing an adiabatic canonical transformation on H\ we obtain the 
transformed Hamiltonian 



a«(t) bf\t) 
c?(t) -a?(t) 



(52) 



where 



fi(t) 



E{ 1} (t) cos 171(f) , b?(t) := friE[ 1} (t) sin 77! (f) 
foEPtysmrntt) , E^it) := £?«(f) cos 77(f) 



ci(f) 



-171 (*) 



^i(f) 

2 / E {1) (s)smr](s)ds . 
Jo 



fi,o :=/i(0)=i/ 



(53) 



Clearly we can repeat this procedure indefinitely and construct an infinite product ex- 
pansion for the evolution operator. Again if we compute only a finite number of terms in 
this expansion, then we obtain a generalization of the adiabatic approximation. The va- 
lidity of this approximation may be checked by computing the transformed Hamiltonians. 
It is not difficult to show that the transformed Hamiltonian obtained after i adiabatic 
canonical transformations is of the form 



Hf\t) = h,{t)S{t) 



where 



h t (t) 
%(f) 



E^\t) cos r](t) cos r]i(t) cos r]2(t) ■■■ cos r]f-i(t) , and 

rt 

2 / E^'(s) cos r](s) cos 771(3) ••• cos ^-i(s) sin r]j_ 1 (s)ds 
J 



where j 6 {2,3, ■■■,£ — 1} and S(t) is a two- by- two matrix of unit determinant. Clearly 
if for some £, he(t) is negligible, then the above mentioned generalization of the adiabatic 
approximation is valid. 

Finally let us note that in general the initial Hamiltonian ([29]) can be written in the 
form 

H{t) = ai(t)a x + a 2 {t)a 2 + a(t)a 3 , (54) 

with oli = (b + c)/2 and a 2 = i(b — c)/2. Performing the canonical transformation (]2lf) 
defined by g(t) = exp{i /q ot2{s)dsa2\ ■, we transform the Hamiltonian ([54]) into 

H'(t) = a'it)^ + a'(t)a 3 = ( °'f } ° ] , where (55) 

\ a (t) — a (t) / 

a'(t) := «i(t) cos^(t) — a(t) sin£(t) , a'(t) := «i(t) sin£(t) + a(t) cos^(t) , 
:= / a 2 (s)^s . 

JO 

Here we have used Eqs. ([H]) and (146|) . Next we perform another canonical transformation, 
namely the one defined by g(t) = exp{i Jq a'(s)dsa3}. This transformation maps the 
Hamiltonian ([55]) into 

e %r] w 
(56) 

where rj'(t) := 2 J * a'(s)ds. This Hamiltonian is not only a member of Class 3 Hamilto- 
nians, but initially (at t = 0) its off-diagonal matrix elements are equal. In particular, 
it has the properties 1. in the above list. Note that we can carry out these canonical 
transformations on any two-level Hamiltonian. Therefore, every two-level Hamiltonian is 
canonically equivalent to a Class 3 Hamiltonian of the form (^). This means that the 
results obtained for Class 3 Hamiltonians apply to arbitrary two-level Hamiltonians. 

V. Time-dependent Simple Harmonic Oscillator 

It is well-known that the solution of every second order linear ODE p3 can be reduced 



to the classical equation of motion for a simple harmonic oscillator with a time- dependent 
frequency uj = u(t), 

x(t) +uj 2 (t)x(t) = . (57) 

It is also well-known that one can reduce both the classical and quantum equations of 
motion for a generalized harmonic oscillator to Eq. (|5"Tp, [23, |25|, EBJ. This equation 



has, therefore, many physical applications p5| , |27| . Yet an exact analytic expression for 
the general solution of this equation is not known even for the case of real frequency |28| . 



B 

In the following we shall consider the case of an ordinary time-dependent harmonic 
oscillator ( 57|) with real frequency. In order to apply the results of the preceding section 



to Eq. (p7|), we first express it in the form of a system of first order ODEs. Defining, 

w*)> : = ( x{t) I > and v (t) ■= A W > 




we can write Eq. (p7|) in the form of the Schrodinger equation ([!]) with a two-level Hamil- 
tonian of the form fl29T) with 



a = , b = i , c = —iuj{t) 2 , and E = uj{t) . (59) 
Since a = 0, this system belongs to the Class 3 of the preceding section with 
/(*) = u(t) , rj(t) = 2 [ u(s)ds, 



= EV{t)( COSr]it l and |f. (60) 

y Uosmr]{t) cos 77 (r) J zujyt) 

Clearly, for real frequency u(t) we can scale the time variable t so that Uq = 1. Then 
the transformed Hamiltonian (|60T) takes the form 



H w (t) = E {1) (t) [sin 77(t)<7! + cos r}(t)a 3 ) = E {1) (fie^®* 3 / 2 ^®**' 2 . (61) 

Note that since rj(t) is real, the Hamiltonian (^1|) is an ant i- Hermit ian matrix with 
orthogonal eigenvectors. Therefore, up to a factor of i it describes a two-level spin system 
with a SU(2) dynamical group [19|, [3|].0 This is rather surprising, for it is well-known that 
the quantum harmonic oscillator has a SU (1, 1) dynamical group and that its Schrodinger 
equation may be reduced to Eq. (|57D by means of a quantum canonical transformation 
corresponding to a time-dependent dilatation |26| . 

In view of the fact that E^ 1 ' is proportional to the derivative of In a;, we can make a 
change of independent variable, namely t — > rj. Note that rj is the integral of a positive 



1 The lack of an exact analytic solution of Eq. (|57]) is not surprising. One way to see this is to recall 
that the time-independent Schrodinger equation for an arbitrary potential V(x) in one dimension is given 

^ + (S^_^)*, = 0, (58) 

where E n and ipn are the energy eigenvalues and eigenfunctions, respectively. Eq. ( |5q ) can be easily 
identified with Eq. ( |57| ) provided that one makes the change of variables: x — * t, ip n — > x, and {h 2 [E — 
V(x)]}/(2m) — > uj 2 (t). This shows that if one was able to find the exact analytic solution of Eq. ( |57| ) for 
arbitrary frequency u>, then one would have been able to find the general solution of the time-independent 
Schrodinger equation for any potential V. 

2 Note that one can absorb the factor of i in the definition of the time variable t, i.e., by defining the 
imaginary time variable t := —it. Therefore, the dynamics given by the Hamiltonian ( |6l| ) may be viewed 
as the dynamics of a spin system with imaginary time. 



real function of t. Hence, it is a monotonically increasing function of t. Making this 
change of variable the Schrodinger equation for the Hamiltonian ( |61~D becomes 

CLT] 

where U(r]) := U'(t(rj)), U'(t) is the evolution operator for the the Hamiltonian (|6T|), 

H(ri) := E(7i)e~ ir,a2/2 cr 3 e iria2/2 = £(77) (sin + cosr]a 3 ) , (62) 

E(rj) := , and J := ^ . (63) 

v/; 2^) 'dry v ; 



Up to a factor of i, the Hamiltonian (|62l) describes the interaction of a spin 1/2 magnetic 
dipole with a changing magnetic field whose direction rotates uniformly in the x — z plane. 

As we mentioned in the preceding section for the Class 3 systems H^ 2 \t) = H(t). 
Hence direct application of the method of the adiabatic product expansion does not lead 
to a solution of the Schrodinger equation for the Hamiltonian (|6l|) or (|62"D. In this case, 
either one constructs the modified adiabatic product expansion of the preceding section 
or examines the adiabatic series expansion of Ref. 0. The latter yields a series expansion 
for the evolution operator U{r]) of the Hamiltonian H(i]), namely 

U( v ) = Te- l tt k{s)ds = l-i P H{s)ds + f P T[H(s 1 )H(s 2 )}ds 1 ds 2 + ■■■ + 

Jo 2 Jo Jo 

(—i) n rv rn 

/ ••• / T[H si (s n )}d Sl ---ds n + --- , (64) 
Jo Jo 



where T stands for the time-ordering operator. Since H(if) is proportional to w'(r]), for 
slowly varying ui one obtains an approximate expression for U{rj) by computing a finite 
number of terms in this series. This is in fact another generalization of the adiabatic 
approximation, because if one keeps only the first term in this series and neglects the 
other terms one is essentially neglecting H or alternatively H^ l > . As we explained above, 
this is just the adiabatic approximation. If one keeps more terms in this series, then one 
obtains a better approximation than the adiabatic approximation. 

VI. Quadrupole Interaction of a Spin 1 Particle with 
a Changing Electric Field 

— * 

Consider a spin 1 particle interacting with a changing electric field £(t) = (£i(t), S 2 {t), £3(1)) 
according to the Stark Hamiltonian 

H{t) = \[j-£(t)} 2 , (65) 

where A is a real coupling constant and J is the angular momentum of the particle. 
The quadrupole interactions of the form (^) have been extensively studied for fermionic 



systems in relation with the non-Abelian geometric phases [29, |3(| |31| (See also [0.) The 
occurrence of non-Abelian geometric phases for the degenerate spin 1 systems has been 
pointed out in Ref. For these systems, the particle has a definite angular momentum 
j = 1 and the Hamiltonian is a 3 x 3 matrix. Using the spin j = I representation of Jj, 
we can express the Stark Hamiltonin (|65D in the form 



H 



' 2 ' 



/ l + 2z 2 y/2. 



ze 



2W 



„2i6 



e 

■\f2ze 
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where r, 9, and z are defined by 



-V2ze ie 



£i + iS 2 



(66) 



2 J 



l + 2z 2 J 



and z := — . 

r 



In view of the general results of Ref. ||33|| , if r ^ then the Hamiltonian (|6q) has a 
degenerate and a nondegenerate eigenvalue. In the following we shall consider the case 
where S3 = 0. The general case £3 7^ can be similarly treated. 

If S3 = 0, then z = and 



H 



v 2 ' 



/ 1 e 
2 

v e 2ie 1 



2/9 



(67) 



The eigenvalues of this Hamiltonian are given by 

E\ = 0, E2 = 



\r 2 . 



(6* 



For r 7^ 0, -Ei is nondegenerate and £3 is doubly degenerate. A set of ortho normal 
eigenvectors of this Hamiltonian is given by 



/ -1 \ 



V e 2ie J 



|lM;J2):=-^ 



/ 1 \ 


„2i0 



\1>2,2;R) ■-- 



1 



(69) 



where i? = (r, 0). 

Next we compute U^ \t) for this system. In order to do this we first use Eq. ( |15| ) to 
calculate A 1 . In view of the fact that the Hamiltonian (67) is Hermitian, \<p n ,a;R) = 
\ip n ,a; R) and Eq. (JT5|) leads to 



A 1 = -d9 , A 2 



-d9 




(70) 



Substituting Eqs. ( [70|) into Eq. (JTJ) and making use of Eq. (j68|) , we find 



^(t) 



3 -i[e(t)-e ] 



K 2 {t) 



-i P (t) 



e -i[e(t)-e ] 
1 



(71) 



(72) 



U&>(t) = \ 



'1 + e~ ip ®)e i6 -® 



(73) 



where 

O ■= 9(0) , and pit) := A / r(s) 2 ds . 

Jo 

Using Eqs. flHf ) and (|69"D, we have 

/ (i + e -ip(t)\ e i0-(t) o (-1 + e~ ip ®)e~®+® \ 

2e-^W o 

v (-1 + e-^wyM*) o 

where 6±(t) :=9{t)±9 . 

Next we compute the Hamiltonian H^'(t). This involves the calculation of A™™(t) 
and H^ n (t) for m ^ n. Using Eqs. ([24] ) and ( |69"D we have 

tr21 



/ 



,21 
^11 



— $ ) ^12 ~~ ^21 ~~ 



e ip(t) ^ #12 = #21 = g 

Substituting these equations in Eq. ( |2"3"D and using Eqs. we obtain 

cosp(t) — isinp(t) \ 



= _0(t) 





i i sin p(t) — cos pit) J 



-9(t)[smp(t)E 2 + cosp(t)S 3 ] , (74) 



where 



J 2 •- 



/ -i\ 


v i y 



and Sq := 



/ 1 \ 


U -l) 



(75) 



It is not difficult to recognize S 2 and S 3 as the Pauli matrices a 2 and 03 represented in a 
(0 + 1/2) representation of SU(2). In view of this identification we can express H^(t) in 
the form 

= -0(t) e ^ El / 2 £ 3 e -^) El / 2 , 



(76) 



where 



Si := 



/ 1 \ 


V 1 o o/ 



and we have used Eq. ([46|). 

The Hamiltonian H^'(t) has the following interesting properties. 

a) For 9 = 0, i.e., 9 = constant, the adiabatic approximation is exact and U(t) = 
U^(t). 

b) In view of Eq. ( |T6"D for 9 ^0, H^(t) has three nondegenerate eigenvalues, namely 
—9, 0, and 9. This is quite remarkable because it shows that the adiabatic canon- 
ical transformation maps the degenerate Hamiltonian ([J7]) into the nondegenerate 
Hamiltonian (ITBT). 



c) Since Sj is a representation of the Pauli matrix <jj, the Hamiltonian flT6|) belongs to 
a representation of the Lie algebra of SU(2). This means that one can reduce the 
Schrodinger equation for this Hamiltonian to that of the dipole Hamiltonian [pi [| 

H dp = -29(t)e lp{t)Jl J 3 e- lp{t)Jl . 

d) We can perform another canonical transformation, namely the one defined by g(t) = 
exp[— zp(i)Ei/2] to transform the Hamiltonian (|76f) into 

= I^Ei - 0(t)E 3 , (77) 

where we have used Eqs. (|76|), ( pT] ) and ([F^). In particular if # and r 2 happen to be 
proportional, i.e., for some cGM 

9{t)=cr{t)\ (78) 

then H^ 1 ' (t) = r(t) 2 (|Si — CE3). In this case the eigenvectors of H^' \t) are con- 
stant and the adiabatic approximation yields the exact solution of the Schrodinger 
equation for (t). The corresponding evolution operator is then given by 

Having obtained the evolution operator for H^' (t) we can use Eq. (^) to obtain 
the evolution operator for H^\t) and H{t). This yields the following expression for 
the evolution operator for H(t): 

U(t) = U {0) (t)e ip{t)Sl/2 U {iy (t) , (80) 

where U {0) (t) and U {1) ' (t) are given by Eqs. (^) and 



The above analysis shows that the condition (|78|) defines a class of exactly solvable time- 
dependent Stark Hamiltonians. If 9 = out, for some constant frequency u, this condition 
corresponds to the case of the rotating electric field E = r (sin ut, cos out, 0) with magnitude 
r. 



VII. Conclusion 

In this article we have extended the method of the adiabatic product expansion to non- 
Hermitian and degenerate Hamiltonians. We showed that in general there were three 
possibilities for the adiabatic product expansion: 

1) The expansion terminates after a finite number of iterations. This happens when 
one of the transformed Hamiltonians vanishes. In this case the method yields the 
exact solution for the Schrodinger equation; 



2) The expansion consists of an infinite number of distinct terms. In this case, the 
method does not lead to an exact solution, but it gives rise to a generalization of 
the adiabatic approximation. This approximation is performed by keeping a finite 
number of terms in the product expansion. The general asymptotic behaviour of the 
adiabatic product expansion has not been studied. However, one can interpret this 
approximation by recalling that the condition for the termination of the product 
expansion corresponds to the validity of the conventional adiabatic approximation 
for one of the transformed Hamiltonians. 

3) The expansion involves terms which are not distinct. In this case the expansion 
does not lead to a solution. However, usually one can make another time-dependent 
canonical transformation after each adiabatic transformation and obtain an infinite 
product expansion with the properties of case 2) above. 

We have considered some specific problems that one can attempt to solve using this 
method. We treated the case of a general nondegenerate two-level system and applied our 
general results to the more specific case of the classical equation of motion for a harmonic 
oscillator with a time-dependent frequency. In this case, we showed that the adiabatic 
canonical transformation mapped the corresponding two-level quantum system to a quan- 
tum system with an anti-Hermitian Hamiltonian. Although the direct application of the 
method of adiabatic product expansion did not yield a solution, we could construct the 
modified adiabatic product expansion. We have also outlined an adiabatic series expan- 
sion for the time-evolution operator of this system which led to another generalization 
of the adiabatic approximation. Finally, we considered the application of our method to 
treat the quadrupole interaction of a spin 1 particle with a changing electric field. The cor- 
responding (Stark) Hamiltonian had a nondegenerate as well as a degerenrate eigenvalue. 
We showed that the adiabatic canonical transformation mapped this Hamiltonian to a 
Hamiltonian which had nondegenerate eigenvalues and belonged to a reducible (0 + 1/2) 
representation of the Lie algebra of SU(2). This means that we can directly use the re- 
sults of Refs. |2], U which treat the Schrodinger equation for a nondegenerate Hamiltonian 
belonging to (an irreducible representation of) the Lie algebra of SU(2). Furthermore, 
we identified a class of exactly solvable spin 1 quadruple Hamiltonians. 
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